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The standard model (SM) invokes the Higgs mecha- 
nism of spontaneous symmetry breaking to generate mass 
for the W and Z bosons, which mediate the weak force. 
The SU(2) x (7(1) symmetry of the electroweak inter- 
action predicts the relation between the W and Z bo- 
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son masses and the electromagnetic and weak gauge cou- 
plings. The prediction for the W boson mass Mw in 
terms of the precisely measured Z boson mass Mz, the 
Fermi decay constant Gf extracted from the muon life- 
time measurement, and the electromagnetic coupling a 
at the scale Mz, is given in the "on-sheh" scheme by [1] 

2 h 3 7ra 1 

W = ~V2G F (l-c 2 w )(l-Ar) ' 

where cyy = Mw/Mz and Ar is the quantum- loop cor- 
rection. A precise measurement of Mw provides a mea- 
surement of Ar. In the SM the contributions to Ar are 
dominated by the top quark and the Higgs boson loops, 
such that Mw in conjunction with the top quark mass 
constrains the mass mn of the undiscovered Higgs boson. 
An mn constraint inconsistent with direct searches can 
indicate the presence of new physics, such as contribu- 
tions to Ar from supersymmetric particles 

The W boson mass ll| has been measured most pre- 
cisely by the LEP ||, |j and Tevatron [5j experiments, 
with the world-average M w = (80392 ± 29) MeV/c 2 [|. 
At the Tevatron, W bosons are mainly produced in quark 
(q') anti-quark (q) annihilation q'q — > W + X. Here X 
includes the QCD radiation that forms the "hadronic re- 
coil" balancing the boson's transverse momentum px @ • 
The W — > tv decays, characterized by a high-p-r charged 
lepton (I = e or fi) and neutrino, can be selected with 
high purity and provide precise mass information. 

This analysis 0i [|[ uses 200 pb _1 collected by the 
CDF II detector 0] in pp collisions at -/i = 1.96 TeV 
at the Tevatron. CDF II is a magnetic spectrometer sur- 
rounded by calorimeters and muon detectors. We use the 
central drift chamber (COT) [|, the central calorime- 
ter [l(| with embedded wire chambers [ll| at the electro- 
magnetic (EM) shower maximum, and the muon detec- 
tors [12| for identification of muons and electrons with 
|ry| < 1 [f| and measurement of their four- momenta. 
The muon (electron) trigger requires a COT track with 
Pt > 18(9) GeV/c [6], and matching muon chamber hits 
(EM calorimeter cluster with Et > 18 GeV). 

In the analysis, we select muons with a COT track 
matched to muon chamber hits and passing quality re- 
quirements, track pt > 30 GeV/c , and a minimum- 
ionization signal in the calorimeter. Cosmic rays are 
rejected using COT hit timing [l3| . We select elec- 
trons with track px > 18 GeV/c , EM cluster Et > 
30 GeV [1, 0] , and passing quality requirements on the 
COT track and the track-cluster matching. Additional 
requirements are based on the ratio of the calorimeter 
energy E to track momentum p (E/pc < 2), the ratio of 
energies detected in the hadronic and EM calorimeters 
Enad/ Eem < 0.1, and the transverse shower profile Q. 
A veto on the presence of a second lepton suppresses 
Z boson background, with negligible loss of W boson 
events. Control samples of Z boson events require two 
oppositely charged leptons with the above criteria. 

The pt of the hadronic recoil (u) equals the vector 
sum u = T,iEiSm(9i)hi/c over calorimeter towers [Io| . 



with energy Ei, polar angle 9i, and transverse direc- 
tions specified by unit vectors hi- Energy associated 
with the charged lepton(s) is not included. We impose 
Pt balance to infer the neutrino's transverse momentum 
p T = | — Pt ~ u\ M an d the W transverse mass mr = 

yfe ( p T p t ~Pt -Pt)/ c - We require p T > 30 GeV/c 
and \u\ < 15 GeV/c to obtain a W candidate sample of 
high purity, whose mj- and lepton pt distributions are 
strongly correlated with Mw- The sample consists of 
63964 W ->■ ev and 51128 W -> fiv candidates. 

The W boson mass is extracted by performing binned 
maximum likelihood fits to the distributions of tut, Pt 
and p T . We generate 800 templates as functions of Mw 
between 80 GcV/c 2 and 81 GeV/c 2 using a custom Monte 
Carlo (MC) simulation [7| of boson production and de- 
cay, and of the detector response to the lepton(s) and 
hadronic recoil. The custom MC optimizes computing 
speed and control of systematic uncertainties. The kine- 
matics of W and Z boson decays are obtained from the 
RESBOS [IH program. We tune the non-perturbative 
form factor in RESBOS, which describes the boson pt 
spectrum at low pt, on the dilepton pt distributions 
in the Z boson data. Single photons (FSR) radiated 
from the final-state leptons are generated according to 
the WGRAD program [lj|. The FSR photon energies are 
increased by 10% (with an absolute uncertainty of 5%) 
to account for additional energy loss due to two-photon 
radiation [l6[ . WGRAD is also used to estimate the initial- 
state QED radiation. We use the CTEQ6M jl7| set of 
parton distribution functions and their uncertainties. 

The custom MC performs a hit-level simulation of 
the lepton track. A fine-grained model of passive ma- 
terial properties is used to calculate ionization and ra- 
diative energy loss and multiple Coulomb scattering. 
Bremsstrahlung photons and conversion electrons are 
generated and propagated to the calorimeter. COT hits 
are generated according to the resolution (« 150 /im) and 
efficiencies measured from muon tracks in T, W, and Z 
boson decays. A helix fit (with optional beam constraint) 
is performed to simulate the reconstructed track. 

The alignment of the COT is performed using a high- 
purity sample of high-p^ cosmic ray muons. Each muon's 
complete trajectory is fit to a single helix [l3|. The fits 
determine the relative locations of the sense wires, includ- 
ing gravitational and electrostatic displacements, with a 
precision of a few microns. We constrain remaining mis- 
alignments, which cause a bias in the track curvature, by 
comparing (E/pc) for electrons and positrons. 

The tracker momentum scale is measured by template- 
fitting the J/ip — > and T — > fifj, mass peaks. The 
J/ip fits are performed in bins of (l/p T {^)) to mea- 
sure any non-linearity due to mismodelling of the ion- 
ization energy loss and other smaller effects, and in 
bins of (cot 0(fi)) to measure the magnetic field non- 
uniformity. To account for the observed momentum 
non-linearity, a downward 6% correction to the pre- 
dicted ionization energy loss is applied in the Simula- 
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FIG. 1: The Z — » /i^i (top) and Z — > ee (bottom) mass fits, 
showing the data (points) and the simulation (histogram). 
The arrows indicate the fitting range. 



tion to make the measured J/ip mass independent of 
(l/p T ((j,)) . The calibration derived from the J/tp and 
T data yields M z = (91184 ± 43 sta t) MeV/c 2 (Fig. Q]) 
from the Z — > jii/z data, consistent with the world aver- 
age [H,|i| of (91188 ± 2) MeV/c 2 . The systematic uncer- 
tainties due to QED radiative corrections and magnetic 
field non-uniformity dominate the total uncertainty of 
0.02% on the combined momentum scale, derived from 
the J/tp, T and Z boson mass fits. 

We simulate the electron cluster by merging the 
energies of the primary electron and the proximate 
bremsstrahlung photons and conversion electrons. The 
distributions of electron and photon energy loss in the 
solenoid coil, and leakage into the hadronic calorimeter, 
are determined using GEANT [l8| as a function of Ex- 
The fractional energy resolution is given by the quadra- 
ture sum of a sampling term (13.5%/ \J Et / GeV) and a 
constant term k = (0.89 ± 0.15)% applied to the cluster- 
energy, and an additional constant term k 7 = (8.3±2.2)% 
applied only to the energies of bremsstrahlung photons 
and conversion electrons. The k 7 term contributes « 
1.3% in quadrature to the effective constant term for the 
inclusive electron sample. The distribution of the under- 
lying event energy [7[ in the cluster is simulated. We tune 
k on the width of the E/pc peak (Fig. [2]) of the W boson 
sample, and k 7 on the width of the Z — > ee mass peak 
when both electrons are radiative (E/pc > 1.06). 

Given the tracker momentum calibration, we fit the 
E /pc peak in bins of electron Et to determine the elec- 
tron energy scale and non-linearity. The position of 
the E/pc peak is sensitive to the number of radiation 
lengths Xo (wl9%), due to bremsstrahlung upstream of 
the COT. We constrain Xq by comparing the fraction of 



E/pc(W->ev) 

FIG. 2: The distribution of E/pc for the W -> ev data 
(points) and the best-fit simulation (histogram) including the 
small jet background (shaded). The arrows indicate the fit- 
ting range used for the electron energy calibration. The jet 
background, which is barely visible on this scale, contributes 
a negligible uncertainty in the calibrations of the electron en- 
ergy scale and the amount of radiative material. 



electrons with high E/pc between data and simulation. 
Applying the E / pc-based energy calibration, we fit the 
Z — > ee mass peak and measure Mz = (91190 ± 67s t at) 
MeV/c 2 (Fig.[T]), consistent with the world average 
For maximum precision, the energy scales from the W 
E/pc fit and the Z — > ee mass fit are combined using 
the Best Linear Unbiased Estimate (BLUE) method [lj| , 
with a resulting uncertainty that is mostly statistical. 

The recoil u excludes towers in which the lcpton(s) de- 
posit energy. The underlying event energy in these towers 
is measured from the nearby towers in W boson data, in- 
cluding its dependence on rj and u. The resolution of u 
has jet-like and underlying event components, with the 
latter modelled using data triggered on inelastic pp inter- 
actions. The recoil parameterizations are tuned on the 
mean and r.m.s. of the pr imbalance between the dilep- 
ton pt and u in Z —> li events. The lepton identification 
efficiency is measured as a function of itii — u ■ p^ /Pt 




« y (W^ev) (GeV/c) | U | (w^iiv) (GeV/c) 

FIG. 3: Left: The M|| distribution for the electron channel 
data (points) and simulation (histogram). Right: The |m| 
distribution for the muon channel. The mean and r.m.s. of 
the histograms agree between data and simulation, within the 
statistical precisions of ~ 1%. 
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using the Z — > 11 data, in order to model its effect on the 
p T and p T distributions. Cross-checks of the recoil model 
using the W boson data show good agreement (Fig.[3|). 
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mr(e,i/) 


80493 ± 48 stat 
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80451 ± 58 stat 
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Pr(e) 


80473 ± 57 atat 
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80349 ± 54 stat 


dz 2Tsyst 


59/48 


PrO) 


80321 ± 66 stat 


± 40 syst 
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Pt(p) 


80396 ± 66 stat 
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44/62 



TABLE I: Fit results and uncertainties for Mw- The fit win- 
dows are 65-90 GeV/c 2 for the m T fit and 32-48 GeV/c for 
the p T and Pr fits. The x 2 °f the fit is computed using the 
expected statistical errors on the data points. 



Backgrounds in the W boson candidate samples arise 
from misidentified jets containing high-pr tracks and EM 
clusters, Z — > 11 where a lepton is not reconstructed and 
mimics a neutrino, W — > tl>, ir/K decays in flight (D1F), 
and cosmic rays. Jet, DIF, and cosmic ray backgrounds 
are estimated from the data to be less than 0.5% com- 
bined. The W — ► tv background is 0.9% (0.9%), and 
the Z -> 11 background is 6.6% (0.24%) in the muon 
(electron) channel, as estimated using a detailed GEANT- 
based detector simulation. The background shapes are 
obtained using simulation and data-derived distributions. 
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FIG. 4: The mr distribution of the data (points) and 
the best-fit simulation template (histogram) including back- 
grounds (shaded), for muons (top) and electrons (bottom). 
The arrows indicate the fitting range. The xVdof for the elec- 
tron channel distribution receives large contributions from a 
few bins near 65 GeV/c 2 , which do not bias the mass fit. 
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39 
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TABLE II: Systematic and total uncertainties in MeV/c 2 for 
the mr fits, which are the most precise. The last column 
shows the correlated uncertainties. 

Table U shows the fit results from the my (Fig. [4j), p T , 
and p T distributions. These fits are partially uncorre- 
cted and have different systematic uncertainties, thus 
providing an important cross-check. The fit values were 
hidden during analysis by adding an unknown offset in 
the range [-100, 100] MeV/c 2 . The systematic uncertain- 
ties (Table |TTJ) were evaluated by fitting MC events to 
propagate the analysis parameter uncertainties to M\y- 



The consistency of the fit results (Table Q]) obtained 
from the different distributions shows that the W bo- 
son production, decay, and the hadronic recoil are well- 
modeled. The statistical correlations (evaluated using 
ensembles of MC events) between the mr and p T (p T ) 
fit values is 69% (68%), and between the p T and p T 
fit values is 27%. We numerically combine (using the 
BLUE method) the six individually fitted Mw values, 
including their correlations, to obtain Mw = (80413 ± 
34 stat ± 34 syst ) MeV/c 2 , with X 2 /dof = 4.8/5. The m T , 
pip and p T fits in the electron (muon) channel contribute 
weights of 32.3% (47.7%), 8.9% (3.4%) and 6.8% (0.9%) 
respectively. This establishes an a-priori procedure to 
incorporate the information from individual fits. The 
muon (electron) channel alone yields Mw = (80352 ±60) 
MeV/c 2 {M w = (80477 ±62) MeV/c 2 ) with x 2 /dof = 
1.4/2 (0.8/2). The m T (p e T , p T ) fit results from the muon 
and electron channels are consistent with a probability of 
7% (18%, 43%), taking into account their correlations. 

In conclusion, we report the first measurement of the 
W boson mass from Run II of the Tevatron. We mea- 
sure M w = (80413 ± 48) MeV/c 2 , the most precise sin- 
gle measurement to date, and we update the world av- 
erage [J to M w = (80398 ± 25) MeV/c 2 . This analysis 
significantly improves in precision over previous Tevatron 
measurements, not only through the increased integrated 
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luminosity but also through improved analysis techniques 
and understanding of systematic uncertainties. As many 
simulation parameters are constrained by data control 
samples, their uncertainties are statistical in nature and 
are expected to be reduced with more data. Inclusion of 
our result in the global electroweak fit 0, 0] reduces the 
predicted mass of the SM Higgs boson by 6 GeV/c 2 and 
decreases its range to ran = 76T. 2 4 GeV/c 2 . 
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